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It is well known that there should be a total cancellation of the IR divergences in unitary inter- 
acting field theories, such as QED and gravity. The cancellation should be at all orders between 
loop and tree level contributions to cross-sections. This is the crucial fact related to the unitarity 
of the evolution operator (S-matrix) of the underlying interacting field theory. In this note we show 
that such a cancellation does not happen in de Sitter space. 
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I. INTRODUCTION 

It is commonly accepted that de Sitter background 
should correspond to the lowest energy state in the grav- 
ity theory with the positive cosmological constant. The 
main argument behind this point of view is that de Sit- 
ter space has the largest possible isometry group with the 
given cosmological constant, while any deviation from de 
Sitter background breaks the isometry. As well it can be 
shown that there are no exponentially growing linearized 
fluctuations in dc Sitter space [1]. 

Let us address, however, the following question in the 
interacting field theory on de Sitter background: Docs an 
inertially moving charged particle in de Sitter space emit 
radiation or not? Because the space in question is confor- 
maly flat, we propose to consider field theory which is not 
conformal, otherwise the behavior of fields is not much 
different from the one in Minkowski space. For example, 
free electromagnetic fields do not feel the expansion of 
de Sitter space and behave as if they are in Minkowski 
one. However, we can consider cither minimally coupled 
scalars and gravitons (whose free field theories already 
are non-conformal) or turn on interactions which break 
conformal invariance. 

It is possible to fetch the answer on the posed question 
even before going into the calculational details using just 
general physical arguments. In fact, an inertially moving 
particle in de Sitter space accelerates with respect to a 
free floating (inertial) observer in the same space. Hence, 
it is tempting to think that inertial particle should radi- 
ate from the point of view of the observer in question. We 
support these general comments with the explicit calcu- 
lation in the main body of the text. Somewhat similar 
phenomenon have been considered in [2]. 

But if the particle radiates it will do that eternally 
— as long as the particle and the background are left 
untouched. Where does the energy for the radiation come 
from? One can object that we should not ask about 
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energy in a time dependent de Sitter like background. 
But exactly in this objection resides the answer to our 
question. In fact, the Hamiltonian of an interacting non- 
conformal field theory in de Sitter background does not 
have a ground state, if the cosmological constant is held 
fixed. We argue that the radiation happens at the cost of 
the decrease of the cosmological constant, if the latter is 
not supported to be constant by any imaginary external 
influence. 

What we are trying to point out here is that the 
situation in many respects is similar to the one in 
QED with the constant electric field if we use the cre- 
ation/annihilation operators which correspond to the ex- 
act harmonics rather than just to plain waves. As is well 
known the similarity goes even further — up to the pair 
creation [3] — but to avoid any mystery of quantum field 
theory in curved backgrounds with horizons, we prefer 
to discuss a simple semiclassical (tree level QFT) phe- 
nomenon and IR rather than UV behavior of quantum 
corrections, which can be completely understood with 
the presently existing level of knowledge. It is obvious 
that, even if we neglect the pair production, a charge 
placed in the constant electric field background will ra- 
diate at the cost of the decrease of the field in question 
or at the cost of work performed by the external source 
keeping the electric field constant. Similarly the cosmo- 
logical constant should decrease if it is not held fixed by 
any external source. 

In this note we would like to find a straightforward 
signal showing that there are problems with quantum 
fields in de Sitter background. Because of asymptotic 
non-flatness of de Sitter space (hence, no energy conser- 
vation) the radiation discussed above could be considered 
as not being a problem, but we show that it inevitably 
results in one. We observe that the evolution operator in 
de Sitter space is not unitary if we keep the cosmological 
constant fixed. The way to see that is through the non- 
cancellation of the IR divergences between tree and loop 
contributions to the cross-sections. 

Let us sketch here the arguments presented in the main 
body of the text. The statement is that in de Sitter 
space a charged particle on mass-shell does emit radi- 
ation. Hence, its virtuality [21] is not related to the 
momentum of the emitted radiation. Thus, the tree 
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level cross-sections for cmmitting soft radiation are fi- 
nite. In fact, recall that IR divergences in cross-sections 
(of a QFT without background fields) appear because the 
propagators of the particles which create radiation, be- 
ing proportional to their inverse virtuality, are singular 
as the momentum of the emitted radiation goes to zero 
[4] : As follows from the energy-momentum conservation 
the virtuality is proportional to the momentum of the 
soft radiated quantum. The power of the singularity is 
such that, after the integration over the phase volume 
of the radiated soft quantum in any cross-section, the 
logarithmic IR divergence appears. Such IR divergences 
cancel with the ones appearing in loop corrections to the 
cross-sections [4] . The cancellation can be directly linked 
to the unitarity of the S-matrix in the QFT [")] or, more 
concretely, to the optical theorem. 

Now in the case of dc Sitter space, while tree level 
cross-sections are finite, the loop diagrams do have IR 
divergences! As a result, unlike the situation in a QFT 
without background fields, there is nothing which can 
cancel them. Hence, the evolution operator in de Sitter 
space is not unitary, as it should be for a non-closed 
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It will become clear from the discussion below that the 
reason for consideration of such a theory is that we would 
like to keep both masses to and M greater than zero. 

This theory in Minkowski space does possess the can- 
cellation of the IR divergences if to = and M > or 
does not have them at all if M, to > 0. As well one can 
consider non-minimally coupled scalars if she will make 
such substitutions as m? rn^+C R for all scalar masses, 
where R is the de Sitter curvature, and C is a parame- 
ter. Conformal coupling corresponds to the case to = 0, 
C = (Z? — 2)/4(Z?— 1). It is important that the interaction 
term breaks the conformal invariance in any case. 

At first sight the most convenient reference frame 
where one can do all the calculations is the planar one 

ds^ — ~dt^ + e^* dxi dxi = — i^—dT^ + dxi dxi) (2) 

where r = e~*, because then we have to deal with the 
non-compact spacial sections and the formulas for am- 
plitudes are very similar to those in Minkowski space 
QED with constant electric field. Unfortunately in this 
coordinates we encounter problems. To see them con- 
sider the Klein-Gordon equation describing propagation 
of free waves in these coordinates: 
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system due to the presence of a background field which 
is held fixed by an external source. 

It is worth pointing out here that the problems with 
cancellation of the IR divergences appear in de Sitter 
space even if we respect the dc Sitter isometry at every 
step of the calculation, i.e. de Sitter space is unstable 
and the isometry is broken, at least if one turns on inter- 
actions. It is unstable in the sense that the cosmologi- 
cal constant will decrease, which will result eventually in 
FRW universe (with non-accelerating expansion) . In the 
latter case we are going to have an analog of the standard 
Minkowski vacuum for quantum fields. 



II. GENERAL DISCUSSION 

We set the cosmological constant to be one and keep 
it fixed throughout the paper. Although our arguments 
are general, for simplicity we would like to consider two 
minimally coupled real scalar particles in D-dimensional 
dc Sitter space with the Yukawa type interaction: 
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and similarly for 'I'. Because this is a free wave equation, 
its solutions obey the superposition principle. Hence, 
from the point of view of the observer, seeing the cor- 
responding metric, particles (single waves) are just so- 
lutions of such an equation, having a finite flux to be 
defined below. But we would like to respect the de Sitter 
isometry, which restricts the choice of the basis of har- 
monics. A particular basis which leads to the de Sitter 
invariant vacuum state is as follows [d]: 

^=y'(^)'^, k^\k\ (4) 

These harmonics correspond to the positive energy 
states, while their complex conjugates — to the nega- 
tive. Here i?^^) jg ^he Hankel function: hI^^ = [ui^^y . 

Let us stress here the main problem with these har- 
monics. The term linear in the differential over r under 
the brackets in (3), which wc refer to as a "friction" term, 
has a wrong sign as t goes to -l-oo (past infinity). As the 
result among the harmonics present in the complete basis 
of the solutions of this equation we have those which are 
exponentially big (in t) in the past, when D > 3. In fact, 
the solution presented in (4) behaves, when r +oo 
(t ^ — cx)), as follows: 
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V'fc (xr^e-*''"^+*'^^ (5) 

This happens because the metric in (2) is singular as r ^ 
+00. As the result all loop diagrams have divergences in 
the r — ^ +00 corner of the time integration axis. Such a 



Here \k,q) = a'^ \vac) and etc.. Here a and b 
are creation operators for the harmonics (4) of the 
fields ^' and ip, correspondingly; \vac) is the de Sit- 
ter invariant Bunch-Davies vacuum [li] and i^i = 

^(^D^f -M^^ = ^{D^-f -m^. Because of 

such a behavior as in (5), the integral in (6) looks like: 

/ + 00 

at the upper integration limit. Hence, for D > 5 we 
have a divergence even in the tree level amplitude inde- 
pendently of the value of m and M. Similar problems 
appear in the loop diagrams starting with _D = 4. 

One can try to avoid the divergence by turning on 
the interactions at some finite tq and evolve to a future 
T < tq. This is explicitly done in loop amplitudes in 
the papers [7] and implicitly — in the papers [s] and ['*]. 
However, in this way we break the de Sitter invariance 
by hand, because the latter acts on tq. Hence, it is not 
an occasion that in the quoted papers a perturbation of 
the de Sitter metric which does not respect the invari- 
ance was observed. It is not that we completely disagree 
with such an approach, taking into account that de Sit- 
ter invariance is going to be dynamically broken anyway, 
but we just would like to show here that one will en- 
counter problems even if she will always try to respect 
the invariance. 

Before going further let us point out the meaning of 
the amplitude (6). First, let us stress that the calculation 
of the amplitude gives a generally covariant (and gauge 
invariant) way to address the question of radiation. In 
fact, if it is not vanishing for given directions of external 
momenta and when all its external legs are on mass- 
shell, it just means that a single wave can split into two 
waves from the point of view of the observer correspond- 
ing to the background metric, in which all calculations 



divergence present along with the IR one which appears 
at r = corner. 

We postpone the discussion of the problems with the 
definition of the S~matrix in de Sitter space to the fol- 
lowing sections and define here the mass-shell three leg 
amplitude as follows. It is proportional to the integral: 
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have been done. 

Second, notice that the amplitude (6) is proportional 
to the spacial ^-function imposing the momentum con- 
servation law p = k + q. In Minkowski space the time in- 
tegral would as well impose the energy conservation law 
Po = ^0 + Qo via (5-function. At the same time on mass- 
shell in Minkowski space we obtain po = ^p^ -|- M^, 

fco = \/P + M2, = \/ cp + rn? and the energy con- 
servation condition does not have a solution. Hence, the 
amplitude is zero and mass-shell (inertial) particles in 
Minkowski space can not emit radiation. 

Now, in (6) we do not have the energy conservation 
due to the presence of the external gravitational field 
originating from the cosmological constant — similarly 
to the QED in an external constant electric field. Thus, 
the amplitude is not zero on mass-shell. For the case 
when it is convergent (i.e. when D < A and M, m > 0) 
this can be seen explicitly via numerical calculation of 
the integral in (6) using Mathematica or Maple. 

In our opinion, the aforementioned divergencies in loop 
diagrams at the past infinity (t = -l-oo) are simply sort 
of boundary effects, which emerged because the planar 
coordinates cover only half of de Sitter space. To avoid 
the aforementioned divergencies in loop diagrams at the 
past infinity (r = +oo) we propose to consider the global 
coordinate system: 



ds^ = -dt'^ + cosh^ tdnjj^i (7) 

where dnl)_^ = ddj + sin^ 6idd'^ + ■■■ + 
SIT? 9i . . .sir? 9£,-2d9\)_^. Unlike the planar coor- 
dinates, the global ones cover de Sitter space completely. 
Important feature of these and the planar coordi- 
nates is that they are seen by inertial observers. The 
Klein-Gordon equation in these coordinates is as follows: 



p ) cx 



oc 



+ 00 



A oc (^k, q 

SD-i)^e^{p-k-<r) (D-i) H^r^ipr) Hl'^Hkr) Hlfiqr) = 
= <5(^-l) [p-k- q) drr^ H^f {pr) Hi]^ {kr) {qr) 
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( + {D-2) t&nhtdt + m2 - ^£^1^] ^ ^ q (8) 
\ cosh t ) 

Here A£i_i(ri) is the Laplacian on the [D — 1)- 
dimensional sphere. 

The "friction" term in (8) is proportional to tanht and 
changes sign from — (D— 1) in the past infinity (t —oo) 
to +(£'—1) in the future infinity {t +oo). As the resuh 
we obtain complete set of modes which are finite at every 
value of t at any given mass m. To solve (8) explicitly 
we can use the separation of variables: 
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where /x = \J rin? — (-^^j^)^ and F[a^b;c] z) is the hy- 
pergeometric function. The solution (10) can be con- 
tinued to the case when m < (Z? — l)/2. The "out" 
modes ffifit) are related to the "in" modes as follows 

The "in" or "out" wave functions in question (as 
well as (4)) are orthonormal with respect to the norm 
i i^i^f—gg^^ do ip2 d^^^fl, which is invariant under 
the change of the spacial section X, as the consequence 
of the equation (8) (or (3)). This norm defines the flux. 
Hence, any solution of the eq. (8) which has a definite 
finite flux (i.e. corresponds to a propagating particle) 
can be decomposed in the complete basis of easer "in" or 
"out" modes. The particular choice of the "in" or "out" 
modes as the basis of harmonics leads to the de Sitter 
invariant vacuum state [11]. 

For the future references let us discuss here the asymp- 
totic behavior of the "in" modes. The hypergeometric 
function F{a, b; c; z) does not have any poles on the neg- 
ative z axis (z = — e^* in our case). Hence, the "in" 
modes (10) are regular at any value of t and m and be- 
have in the past infinity (i — > — c») as 

<^± ^e(^T»A^)t (11) 

because F 1 as 2 = — e^* 0. In the future infinity 

{t +00, z = — e^* —00) they look like 

vf^e-^*{c,e^'>'' + C2e^"'') (12) 

with some complex constants ci and C2. Such a behavior 
follows from 



ijj,,{t,n)^^j{t)Y,n{^) (9) 

Here A^i-i (17) „(17) = -j(j + D - 2) Y,„(17), and n is 
the multi-index (rii, . . . , iio-i)- 

The spherical harmonics Yj „(r2) have obvious proper- 
ties presented in [10]. The field (pj{t) obeys the obvious 
equation following from (8) and (9). This equation has 
two designated complete sets of solutions: so called "in" 
and "out" modes [11] (see as well [10]). The complete set 
of "in" modes is 



(10) 



r 



lim F{a,h;c-z)^ci{-z)-'' + C2{-z)-^ (13) 

z — * — 00 

Note that if m = we have harmonics which approach 
non-vanishing constants as t ^ ±00, but there are no 
any modes which are exponentially growing. 

For the future reference let us define here the propa- 
gator for the "in" modes [11], [10]. Consider de Sitter 
invariant function depending on two points [11], [10]: 

Z(z,z') = — sinht sinht' + coshi cosht' cosA17 (14) 

where Af2 is the angle between the spacial parts of the co- 
ordinates z and z' . It can be shown that Z(z, z') = cos L, 
where L is the geodesic distance between z and z' for 
spacial separations, or i times the geodesic proper time 
difference for time-like separations. The Green function, 
being de Sitter invariant, should depend only on such a 
combination of the two points. Hence, the equation for 
the Green function, which is (8) with the appropriate 
i5-functional sources on the RHS, can be converted into: 

[(1-^2) dl^DZdz~m^] G{Z) = 

= A5{Z + 1) + B5{Z -1) (15) 

by the direct change of variables (14), from {t,Vl) to Z. 
Here A and B are some constants. The RHS of this 
equation is singular when the z and z' points coincide (i.e. 
when Z = 1) and when z coincides with the antipodal 
point of z' (i.e. when Z ~ —1) [11]. The "in" Feynman 
type propagator obeys this equation with A = B = \ 
and is given by [11]: 



^f{t) oc cosh^(0 e(^+-T^T"^) * F (^j + j + 



D - 1 
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G„(Z) cx 



sinh(7r /i) 



F 



F 



D - 1 
D- 1 



D-1 



D 1 



2 

D- 1 



L> 1 - Z 



(16) 



where /i is defined above. We are going to use this func- 
tion in the loop calculations below. 
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III. ON QFT WITH COMPACT SPACIAL 
SECTIONS 



We see that spacial sections in global coordinates are 
compact {D — 1) -dimensional spheres, which is less con- 
venient than flat sections in planar coordinates. To see 
that compactness of spacial sections does not spoil all 
the picture, in this chapter we would like to consider a 
general features of QFT on such a space-time as, for ex- 



ample, ds 



with fixed radius R. We 



are going to show that such a QFT as (1) on the back- 
ground in question has similar properties to the theory 
in Minkowski space: such properties as the impossibility 
for inertial particle to emit radiation as measured by an 
incrtial observer and the cancellation of IR divergences. 
Let us stress here that our confidence in the fact that 
there are no problems with the theory (1) on the back- 
ground in question is relaying on the obvious observation 
that it has the unitary evolution operator. 

Definite energy mass-shell harmonics in such a the- 



ory look like 



oc 



-i ko t 



yjm(ri), where fco 



'm j {j + D — 2). Then, the three leg mass-shell 



amplitude is: 



+ 00 



J — OO 

dnY,,„Am*rjm*n,m 



(17) 



With such pq, ko and qo the condition po = + Qo can 
not be saturated for 72 — js < ji < j2 + js- Hence, 
the argument of the (S-function imposing the energy con- 
servation is always non-zero, i.e. the amplitude itself is 
zero. Thus, similarly to the Minkowski space, we arrive 
at the obvious conclusion that in the space in question 
an inertial (mass-shell) particle can not emit radiation. 

Recall that in Minkowski space IR divergences in the 
cross-section appear only when m = 0. Amplitude of 
a hard process containing emittion of one soft ip mass- 
shell quantum by the heavy particle has the three leg 
multiplicative contribution: 



dt / dilGAiit', n'; t, n) 

XD J 



(18) 



where Gm is the propagator of the field, i.e. one of 
the legs in the amplitude is off-shell. The propagator is 



Git, n; t', n') = J2 



^x{t, n)'fi{t', o') 



(19) 



A 



where [D-AP]'^x = A*a, and □ = 
^Ad-i(^). Obviously *a oc e~''">*Yj.,n2{^), where 
now fco = sJm^ + \+^J2 {j2 + D-2), i.e. = 

J dko X]j2=o- ^^^^ integral over t in (18) leads to the 
energy conservation of the form: 



and is proportional to S{po — fco ^ Qo)- The second in- 
tegral (over the angles) gives the generalization of the 
3j symbols to the SO{D — 1) group with D > A, which 
is not quite convenient object in comparison with the 5- 
function appearing for the case of the flat spacial sections. 
The 3j symbols are nonzero if j2 — J3 ^ ji ^ J2 + J3 ■ 
On mass-shell we have that 



D-2), 



A/2 



i?2 



J2 ij2 +D-2), 



^^' + ^+]^^2 {j2 + D-2) 



+ Us + d-2) 



(20) 



For the big ji ~ j2 and small (soft) we have the solu- 
tion of this equation as follows: A 

where po — 



M2 



-2 (po g o - ji J3) 

ji {ji +D-2) and 
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^ js (-D — 2). Hence, the amplitude is divergent as 



1/A cx l/\/j3 when 0, while the cross-section is 
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divergent as l/js- Similar divergences (with the oppo- 
site sign) appear in one loop contributions to the cross- 
sections. 

To clarify the situation let us consider the Minkowski 
space (ds^ = —dt^+dx^) variant of the theory in question 
and the part of an amplitude responsible for the radiation 
of a soft tp particle by the hard "if one. For simplicity here 
we restrict ourselves to the four space-time dimensions. 
If we insert the Minkowski space analog of the propaga- 
tor (19) into the amplitude, we can use the (5-functions, 
imposing energy-momentum conservation at the vertex, 

to fix A = fc2-fc2_^2 ^ (pp_^jj)2_(p._-)2_^2^ 

is the four-momentum of on-shell incoming particle, q 
is the four-momentum of the on-shell outgoing radiated 
■0 particle and k is the four-momentum of the off-shell 
(with virtuality A) outgoing, after the radiation, ^' par- 
ticle. 

If we consider radiation of the very soft particle, i.e. 
the modulus of the corresponding D-momentum is \q\ 
0, then k is very close to the mass-shell, while = M"^ 
and = 0, i.e. A = —2pq, because m = 0. Thus, 
the propagator is singular as \q\ —^ 0. Moreover, such a 
dependence of A on g is important for the factorization of 
the IR divergences in the cross-sections for radiation of 
many soft quanta [12], which, in its own right, is crucial 
for the total cancellation of all divergences. 

As the result, after the integration of the differential 
cross-section for the radiation of one soft quantum over 
its invariant phase volume, we obtain: 

This is the IR divergence with the cutoff mo 0. Be- 
cause all IR divergences are of the same order, we have to 
sum such contributions over all external legs of the hard 
process in question [4], [12]. Then similar divergences 
(with the opposite sign) appear in the loop contributions 
to the cross-section of the hard process: 

loop IR divergence cx / d^q- — -5 — ^. 

J (PQ) q 

All such contributions (from loops and tree level dia- 
grams) add up, so that every divergence does cancel [12]. 
Higher loop contributions cancel with the divergences 
coming from multiple soft quantum radiations [4]. It 
is important to stress here that we can choose another 
basis of harmonics for 5* (dressed with the cloud of ip's) 
such that IR divergences will not be present neither in 



trees nor in loops [5], but it is impossible to get rid of 
the divergences only in trees without cancelling them in 
loops or vise versa. 

IV. RADIATION AND IR DIVERGENCES IN 
DE SITTER SPACE 

In this section we are going to show that in de Sit- 
ter space the IR divergences do not cancel already at 
the leading order. But before going into the calcu- 
lational details let us note that we avoid using the 
term S-matrix in de Sitter space. The latter is build 
on the basis of matrix elements of the evolution op- 
erator S = T e^' ^-^ '^^ where iJ,„t is the in- 
teraction Hamiltonian. The matrix elements in ques- 
tion are defined with respect to the basis of states 
{out \ a ... a and a'*' . . . \in) while it is assumed that 

\out) = e~* -^-^^ \in) = (phase) \in), where Hq is the 
free Hamiltonian. In the odd space-time dimensional de 
Sitter spaces \out) = \in). Hence, in odd dimensions we 
are safe and our arguments pass smoothly [lo]. 

But it appears that in even dimensional de Sit- 
ter spaces \in) ^ \out) [11], [10]. In fact, 

in even dimensions \out) = e^^ ■^-'^ ^° \in) =/= 
(phase) I m), because \in) is not an eigen-state of the 
free Hamiltonian. The latter has the form Hq oc 
J2k [Mt) flfc «fc + B{t) d-k dk + B*{t) dtk flfc ] with some 
functions of time A{t) and B{t) [l.'i]. Exactly due to the 
presence of the non-diagonal terms d-k dk and dl^ d^ in 
the Hamiltonian we have to make the Bogolyubov trans- 
formation to diagonalize it and to observe the particle 
production [3]. 

Thus, \out) state differs from \in) by the presence of 
the created particles, which can be explicitly established 
by the following relation \out) = V^(d,d+)|m) with some 
operator V [13]. In this note we would like to consider 
the matrix elements of the evolution operator which are 
of the form \ d . . .dS d^ . . .d^ \in) . Physically this 
means that wc neglect the particle production by the 
external field and consider only scattering amplitudes in 
such a background. Wc strongly believe that this is suf- 
ficient to make a statement about (non~)unitarity of the 
evolution operator S itself. To see that our arguments 
are meaningful one can consider the similar situation ap- 
pearing in QED in the background of the constant electric 
field. 

It is straightforward to find the basic tree level mass- 
shell amplitude describing the process when \E' particle 
radiates the particle. In the global coordinates the 
amplitude is proportional to: 
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J J —oo 



D - 1 



2 

L> - 1 



(21) 



where /ii = - (i^)^ = - 

Note that (21) is vahd even if m or Ad are less than 
{D — l)/2. Such an ampUtudc is just an analytical con- 
tinuation of the corresponding generalized 3j~symbol, 
which follows from the continuation 5^ — *■ dS^ . From 
this we can already argue that the mass-shell amplitude 
(21) is non-zero. 

However, if either one of the masses, m or M , is van- 
ishing, the integral for the amplitude is divergent (see 
below). We discuss the meaning of these divergences in 
the concluding section. At this stage we would like to 
avoid such problems with divergences of the tree level 
amplitudes and to keep our discussion as transparent as 
it is possible for the case in question. It happens that, 
if we keep both masses M and m non-zero, the integral 
(21) is convergent. Unfortunately, even Mathematica and 
Maple refuse to take such an integral analytically. Let us 
show explicitly that it is really convergent. 

The integrand expression in (21) can hypothctically 
blowup only if t ±00, because the hypergeometric 
function is regular for the negative argument, i.e. for any 
finite value of t. As < — > —00 we can use the behavior of 
the "in" harmonics from (11) to obtain that the integrand 
expression in (21) approaches «A'2)t fgj- ^j^g lower 

integration limit. Hence, the integral is convergent in 
this corner of the integration axis even if is purely 
imaginary, i.e. when m < {D — l)/2. Indeed in the latter 
case |/i2| < (-D — l)/2. This inequality is saturated only 
when m = 0. It is only in the latter case there can be 
the perfect cancellation of the exponential suppression 
e'^~ «M2)t _> X as i ^ —00, and we have the divergent 
amplitude. 

In the other corner of the integration axis, i.e. when 
t ^ -I- 00, we can use the asymptotics as in (12) and find 
that the integrand behaves as 



■t (cie+*''i* 



-i /^i t 



Co e 



C2 e 



(22) 



Hence, in this corner of the t integration axis the integral 
(21) is convergent if fii is real and m ^ 0. 

As well, using Mathematica and Maple for numerical 
calculation of the integral (21), one can explicitly sec that 
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FIG. 1; Real part of the integrand in (21) for different values 
of M, m and ji, j2, ja- Imaginary part is equal to zero for 
the first three plots. The presence of non-zero imaginary con- 
tribution to the leading amplitude is already a sign favoring 
that the theory in question is non-unitary. 



it is not zero for j2 — js < ji < j2 + ja ■ The integrand 
expression in (21) is plotted on Fig. 1 for several different 
values of A/, m and ji, j2, js- Thus, all our considerations 
so far at least mean that a massive particle can radiate 
another massive particle on mass-shell in de Sitter space. 

Based on these considerations we can make a general 
conclusion that mass-shell particles can radiate fields 
under which they carry charges, unless the correspond- 
ing theory, describing interactions between "matter" and 
"radiation" , is conformal. For example, all particles can 
radiate gravitons, we just have to appropriately under- 
stand the corresponding divergent amplitudes as the gen- 
eralized functions. With the similar reasoning one can ar- 
rive at the conclusion that eternally accelerating charged 
particle in Minkowski space (e.g. under the action of the 
constant electric field) does emit radiation. 

As the side remark let us point out here that one may 
object the conclusions made in the previous paragraph 
based on the following considerations. It is known that if 
one drops a spherically symmetric massive (m) body in 
de Sitter space it produces the static de Sitter black hole 
metric outside its own volume: 



8 



ds^ = - 1 - — dt^ + 7 ^ ^ 

+r^dnl,_2 (23) 

and docs not produce any (non-static) gravitational 
waves on top of that. Hence, it seems that this argu- 
ment precludes our conclusion that an inertial massive 
body in de Sitter will produce gravitational waves. But 
the important point is that such a metric as (23) is seen 
by a non-inertial observer which is fixed above the sur- 
face of the spherical body, i.e. the body and the observer 
compose a bound state and do not move with respect to 
each other. Thus, it is not an occasion that such an ob- 
server does not see any radiation from the massive body. 
At the same time our statement is that it is the iner- 
tial observer which sees the radiation from free floating 
bodies in de Sitter space. 

Let us see now what happens with the cross-section 
of a hard process containing the radiation of the soft ip 
quantum by in de Sitter space. Because the amplitude 
(21) is non-vanishing on-shell the virtuality of the ^' 
particle is not related to the mass-shell momentum (ja) 



of the radiated V' particle. Indeed, to obtain the multi- 
plicative factor in the amplitude of a hard process, corre- 
sponding to the radiation of the soft quantum by the hard 
one, we have to multiply by 1/A the same amplitude as 
(21) with the only change of in the second wave func- 
tion under the integral in (21) by yM2~rA~~(^^i)^. 
Obviously the amplitude is not singular as 0, be- 

cause A is not related to j^. In fact, we do not have S- 
function (imposing energy conservation) which fixes the 
value of A as it was in Minkowski space. Hence, we just 
have to integrate over all possible values of the virtuality 
in the amplitude, unlike the Minkowski space case. Thus, 
the cross-section is not divergent as well, which can be 
explicitly seen via similar reasoning to the one presented 
after cq. (21). But even if it was divergent, we would 
not have had the factorization of the divergences due to 
such a behavior of A. The latter fact anyway spoils the 
cancellation of the IR divergences at all orders [12]. 

At the same time it happens that loop diagrams in 
de Sitter space do have IR divergences even for massive 
particles. Consider the one loop self-energy diagram for 
the ^ particle. It has the contributions of the form: 



/+00 r+oo p p 

dti / dt2 / dfli / d02cosh'°"^(ti) cosh^~^(t2) x 

D-l D-1 D 1±Z\ 



xF 



D - 1 



D 1±Z 



1 112- — ; 

' ' 2 2 



(24) 



where Z is given by (14) and we have borrowed prop- 
agators from (16). It is straightforward to see (using 
(13)) that such an integral divergent in the IR, i.e. as 
Z(z, z') -^■ oo (see e.g. similar discussion in [14]). In- 
deed, if say ti — > — oo, while t2 +oo, then according 
to (14), Z{z,z') e*2~*i (H-cosrj)/4 and the integrand 
expressions in (24) behave as: 



2[ e~*''^ (t2-tl) _|_ ^'^ giM2 (*2-tl) 



(25) 



Hence, the integrals over ti, t2 in (24) are divergent if 
fi2 is pure imaginary, i.e. when m < (D — l)/2. Fur- 
thermore, we have the IR divergence in the causally con- 
nected region, because any points with t —^ —oo and 
t +00 are causally connected in de Sitter space. 
Hence, restricting oneself to the region within the cos- 
mological horizon does not help to cut or get rid of such 
IR divergences. 



Thus, for any Af > 0, but < m < (£) - l)/2, we 
obtain finite tree level contributions to cross-sections and 
there is nothing which can cancel the IR divergences in 
the loop diagrams. It is probably worth pointing out here 
that we can not interpret the divergence in (24) as an 
analog of the coUinear one [4] for many obvious reasons. 
At least it is present for any value of the mass M of the 
particle emitting the radiation. 



All the considerations above make us to conclude that 
the evolution operator leading to such a diagram technic 
is not unitary and the system of de Sitter background 
plus QFT is not closed if the cosmological constant is held 
fixed. It is important to stress that in the circumstances 
under consideration one can not make the IR divergences 
to cancel by a unitary change of the basis of the creation- 
annihilation operators. 
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V. CONCLUSIONS 

Before drawing any conclusions let us make a few side 
remarks about the divergences of (21) when either of the 
masses M or to is vanishing. Similar divergences ap- 
pear in anti-de-Sitter space: note the similarity of the 
anti-de-Sitter metric ds^ = ^ {dz^ + dxa dx'^) to the 
one in (2) with the crucial difference, however, that z 
is not time-like. Hence, anti-de-Sitter space is not glob- 
ally hyperbolic (because there is time-like boundary) but 
does not have an event horizon (because there is a glob- 
ally defined time-like Killing vector) . Because of the lack 
of the global hyperbolicity, which results in such a well 
known effect that in anti-de-Sitter space waves can repel 
from the spacial boundary, one can not define Cauchy 
problem in such a space [15]. Thus, while it is possi- 
ble to define the unique anti-de-Sitter invariant vacuum 
state, one can not define appropriately the evolution op- 
erator for a QFT in such a background. According to 
the AdS/CFT-correspondence [16], [17], [18] one treats 
the IR divergences in the QFT on anti-de-Sitter space, 
which appear in the wave functional rather than in the S- 
matrix, as the UV divergences in the QFT on its bound- 
ary. 

On the contrary, although we do not have Poincare 
invariance in de Sitter space, it is globally hyperbolic and 
one can define evolution operator there. Let us stress here 
that we do not have Poincare invariance, for example, in 
the presence of the constant electric field in QED, but 
we still can define the evolution operator, because the 
evolution problem can be correctly formulated in such 
circumstances. 

At this point we can and should address the question 
why anti-de-Sitter space is stable? It seems that a free 
floating particle in anti-de-Sitter space will emit radi- 
ation as well. However, this question and the question 
of the cancellation of the IR divergences can not be for- 
mulated in anti-de-Sitter space because of the impos- 
sibility to define the time evolution operator in such a 
background due to the lack of the global hyperbolicity. 

Let us now come back to the conclusions. We see that 
QFT in de Sitter space (i.e. if we fix the cosmological con- 
stant) behaves as if it is formulated in a background of 
an external quasi-classical gravitational field (excitation 
above the correct vacuum) due to cosmological constant 
— analogously to the QED in a constant (in space and 
time) electric field, i.e. as the non-closed system. That 
is the only interpretation of the non-unitarity, which we 
can give. Hence, the only conclusion which we can make 
here, without performing a direct calculation, is that the 
cosmological constant will relax to zero creating particles 
via Gibbons-Hawking pair production (which we do not 
discuss in our paper) and via performing work on accel- 
erating created particles, which leads to the radiation in 
its own right. 



The conceptual question is how fast will be the relax- 
ation of the cosmological constant? If the relaxation goes 
fast enough we can hope to explain via this mechanism 
the cosmological constant problem along with obtaining 
natural inflation without any inflaton field [7]. We think 
that the rate of the decay should just depend on the ac- 
tual magnitude of the field. For the big enough field the 
rate of the decay should be big: bigger the energy pool is 
— easier to create light particles. The technical question 
is how to calculate the rate of the decay of the cosmolog- 
ical constant? 

If the cosmological constant is very big we have to use 
theory of quantum gravity. Unfortunately string theory 
does not have a formulation on de Sitter space, because 
such a background spoils conformal invariance of the the- 
ory on the string worlds-sheet. In this case we have to 
deal with off-shell formulation of the closed string the- 
ory, which is very hard to do with the presently existing 
first quantized variant of the theory (see, however, [14]). 

Thus, the only hope is that we can use the ordinary 
Einstein-Hilbert theory, which dominates in the IR, if the 
cosmological constant is smaller than the Plank scale. In 
this respect we should stress that the issue of the insta- 
bility of de Sitter space, of the IR divergences and of the 
non-unitarity have been discussed in various places [2], 
[7], [11], [14], [19], [20]. As well the running of the cos- 
mological constant due to the quantum fluctuations have 
been found in [S], [7] and [!)] using either Heisenbcrg or 
Schwinger-Keldish technics. Here, apart from presenting 
the listed above new phenomena supporting the conclu- 
sion that de Sitter space is unstable, we should criticize 
the actual calculation of the decay rate of the cosmolog- 
ical constant performed in [8], [7] and [9], because, as 
we just pointed out, the rate was found with the use of 
the non-unitary evolution operator. Hence, the questions 
posed above remain to be answered. 
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